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Abstract 

A vibrant and healthy gut flora is essential for preventing the proliferation of Clostridium difficile, a pathogenic bacterium 
that causes severe gastrointestinal symptoms. In fact, most C. difficile infections (CDIs) occur after broad-spectrum antibiotic 
treatment, which, by eradicating the commensal gut bacteria, allows its spores to proliferate. Hence, a C. difficile specific 
antibiotic that spares the gut flora would be highly beneficial in treating CDI. Towards this goal, we set out to discover small 
molecule inhibitors of the C. difficile enzyme dehydroquinate dehydratase (DHQD). DHQD is the 3 rd of seven enzymes that 
compose the shikimate pathway, a metabolic pathway absent in humans, and is present in bacteria as two phylogenetically 
and mechanistically distinct types. Using a high-throughput screen we identified three compounds that inhibited the type I 
C. difficile DHQD but not the type II DHQD from Bacteroides thetaiotaomicron, a highly represented commensal gut bacterial 
species. Kinetic analysis revealed that the compounds inhibit the C. difficile enzyme with K, values ranging from 10 to 20 uM. 
Unexpectedly, kinetic and biophysical studies demonstrate that inhibitors also exhibit selectivity between type I DHQDs, 
inhibiting the C. difficile but not the highly homologous Salmonella enterica DHQD. Therefore, the three identified 
compounds seem to be promising lead compounds for the development of C difficile specific antibiotics. 
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Introduction 

The human pathogen Clostridium difficile is a spore-forming, 
Gram-positive, anaerobic bacillus that secretes two types of toxins, 
which induce severe diarrhea, fever, and nausea. Notably, the 
number one risk factor for C. difficile infection (CDI) is treatment 
with a broad-spectrum antibiotic to combat a preexisting bacterial 
infection [1]. In fact, 15-25% of all antibiotic-associated diarrhea 
cases are caused by C. difficile [2,3]. Under normal conditions, the 
myriad of bacterial species in the gut successfully outcompete C. 
difficile. However, once broad-spectrum antibiotics are adminis- 
tered, the gut flora is largely eliminated, allowing C. difficile spores 
to proliferate without competition from the non-pathogenic 
bacteria [1]. In other words, the normal gut flora acts as an 
anti-C. difficile environment. Hence, physicians treating patients 
with CDI are challenged by two opposing goals. One goal is the 
eradication of C. difficile and the original bacterial infection for 
which antibiotic treatment was initiated. This would be achieved 
by continued administration of broad-spectrum antibiotics, with 
the downside of concomitant destruction of the gut flora. The 
other goal is cessation of antibiotic treatment to allow the recovery 
of the gut flora that is required to combat C. difficile, with the 



downside that ending antibiotic administration may allow C. 
difficile to proliferate in the period prior to flora recovery. This 
Catch-22 scenario could be resolved with a C. difficile specific 
antibiotic, which would prevent C. difficile proliferation while 
allowing for the repopulation of the gut by commensal bacteria. 
This concept is validated by fidaxomicin (Dificid), the first C. 
difficile-targeting narrow-spectrum antibiotic recendy released to 
market, which has been shown to reduce the likelihood of CDI 
reoccurrence [4] . However, fidaxomicin affects not only C. difficile, 
being active against many Gram-positive bacteria, especially 
Clostridia [5] . This fact, and the likely emergence of fidaxomicin- 
resistant C. difficile, indicate a clear need for additional highly C. 
difficile specific antibiotics. 

Towards the goal of developing a narrow-spectrum agent for 
CDI, we commenced studies of shikimate pathways enzymes. The 
ultimate product of this 7-step pathway is chorismate, a precursor 
required for the biosynthesis of the three aromatic amino acids as 
well as other important metabolites. Because humans lack the 
pathway and must obtain the aromatic amino acids through 
dietary sources, the enzymes involved in shikimate biosynthesis 
provide suitable targets for antibacterial drug discovery [6]. The 
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3 r -step of the shikimate pathway involves the conversion of 3- 
dehydroquinate (DHQ) to 3-dehydroshikimate (DHS). Interest- 
ingly, the enzymes that catalyze this reaction, dehydroquinate 
dehydratases (DHQDs), are represented in bacteria by two 
different subtypes, I and II [7,8]. Present in C. difficile, type I 
DHQDs are ~30 kDa enzymes that assemble into homodimers 
[7] and use a Schiff-base intermediate to catalyze the dehydration 
reaction [9] . Present in the highly prevalent commensal bacterial 
species Bacteroides thetaiotaomicron and Bifidobacterium longum, among 
others, type II DHQDs are ~ 1 7 kDa in size, form homododeca- 
mers, and catalyze the reaction via an enolate intermediate [7,10]. 
Since type I and type II DHQDs have different binding sites and 
reaction mechanisms, we reasoned that it should be possible to 
identify small molecules that selectively inhibit the type I C. difficile 
DHQD (cdDHQD) but not the type II B. thetaiotaomicron DHQD 
(M3HQD) and that such an inhibitor would have the desirable 
property of inhibiting C. difficile proliferation while being 
compatible with continued growth of a large subset of the 
commensal bacteria. 

Here we present the discovery and characterization of three 
caDHQD inhibitors identified by HTS. Surprisingly, in addition to 
exhibiting anticipated inter-type selectivity, the inhibitors also 
exhibit intra-type selectivity, inhibiting the cdDHQD but not the 
homologous Salmonella enterica type I DHQD (.veDHQD). This work 
provides a proof-of-principle, demonstrating the feasibility of 
selectively targeting «DHQD while sparing homologous enzymes 
from non-pathogenic bacterial species and opens the path for the 
development of a novel class of C. difficile specific antibiotics. 

Materials and Methods 

Gene Cloning and Enzyme Expression and Purification 

Clostridium difficile aroD (cdDHQD), Salmonella enterica 
aroD (seDHQD), Vibrio cholerae aroE (V. cholerae SDH), 
Bacteroides thetaiotaomicron aroK (B. thetaiotaomicron SK), 
Bacteroides thetaiotaomicron aroQ (btDHQD), Vibrio cholerae 
aroQ (vcDHQD), and Yersinia pestis aroQ (ypDHQD) were 
amplified from genomic DNA by PCR and subcloned into the 
pMCSG7 expression vector. The BL21 (DE3) E. coli strain was 
used for recombinant expression for all but btDHQD (type II 
DHQD), which was expressed in the KRX E. coli strain since this 
enzyme was insoluble in BL21 cells. For expression, 1-3 liters of 
TB media were inoculated with appropriate starter culture for 
each protein and shaken at 225 RPM at 37°C. When an optical 
density of ~0.8 at 600 nm was achieved, protein over-expression 
was induced by adding isopropyl- 1 -thio-D-galactopyranoside to a 
concentration of 0.5 mM, the temperature was reduced to 25°C, 
and the culture was left overnight. The following morning, cells 
were harvested by centrifugation and lysed by sonication in a 
buffer containing 10 mM Tris (pH 8.3), 500 mM NaCl, 10% 
glycerol, and 5 mM P-mercaptoefhanol. The resulting lysate was 
cleared by centrifugation, loaded onto a 5 mL His-Trap HP Ni 
Sepharose column (GE Healthcare), washed with a buffer 
containing 10 mM Tris (pH 8.3), 500 mM NaCl, 25 mM 
imidazole, and 5 mM P-mercaptoethanol, and eluted in a buffer 
containing 10 mM Tris (pH 8.3), 500 mM NaCl, 500 mM 
imidazole, and 5 mM (3-mercaptoethanol. The resulting elutant 
was injected onto a S-200 gel filtration column (GE Healthcare) 
equilibrated with buffer containing 10 mM Tris (pH 8.3), 
500 mM NaCl, and 5 mM (5 -mercaptoethanol. For each 
purification, SDS-PAGE chromatography confirmed that the 
major peak off the gel filtration contained a single major band 
consistent in molecular weight with that anticipated for the 
recombinant protein. To remove expression tag, cdDHQD, 



seDHQD, and btDHQD were incubated overnight at 4°C with 
hexa-histidine tagged TEV protease, reloaded onto the 5 mL His- 
Trap HP Ni Sepharose column, and eluted with 25 mM 
imidazole. 

Cc/DHQD High-throughput Screen 

Automated high-throughput screening against a library of 
50,000 Chembridge compounds was performed at room temper- 
ature on a Tecan Freedom EVO 200 liquid handling platform 
with an integrated Tecan Infinite F200 Pro microplate reader. All 
assays were performed in duplicate in black, flat-bottom 384-well 
plates containing a final reaction volume of 50 |jL. The assays 
were assembled as follows: 40 |0,L of enzyme solution (50 mM 
HEPES, pH 7.5, 0. 1 mg/mL BSA, 5 mM DTT, 0.01 % Triton-X, 
1-5 nM «/DHQD, 125 nM V. cholerae SDH, 125 nM B. 
thetaiotamicron SK) was dispensed into wells and incubated with 
200 nL of 10 mM inhibitor (40 pivl final concentration) for 5 min. 
Reactions were then initiated with 10 |iL of substrate solution 
(50 mM HEPES, pH 7.5, 0.1 mg/mL BSA, 5 mM DTT, 0.01% 
Triton-X, 25 mM MgCl 2 , 5 mM ATP, 125 uJVI NADPH, 
250 |jM DHQ, shaken vigorously for 20 s, and measured for 
fluorescence emission intensity (excitation X: 340 nm, emission X: 
460 nm) over a period of 7 minutes. Slopes from t = 1-5 min were 
recorded for each well. Each 384-well plate contained 32 positive 
control wells (200 nL of DMSO replacing 200 nL of compound in 
DMSO) and 32 negative control wells (assay components lacking 
cdDHQD). Compounds producing &40% inhibition were selected 
for further analysis. 

Cc/DHQD Hit Validation and Selectivity Experiments 

Putative hits identified in the primary screen were cherry- 
picked, retested in the HTS primary assay, and counter-screened 
against the SDH/SK coupling enzyme reaction. The counter- 
screen was assembled in the same manner as the primary assay, 
with the following changes: cdDHQD enzyme was eliminated from 
the enzyme solution, the concentration of SDH was decreased to 
1 nM, and the DHQ substrate was replaced with DHS for a final 
concentration of 50 |0M. Compounds that reproducibly inhibited 
cdDHQD and that displayed no inhibition against the coupling 
enzymes SK and SDH at 40 |jM compound were then assayed 
against .seDHQD and feDHQD using the primary screen DHQD 
assay and up to 200 |0,M compound. The three confirmed hits 
were repurchased from Chembridge and validated against 
cdDHQD in an alternative assay monitoring the conversion of 
DHQ to DHS at 234 nm (50 mM potassium phosphate buffer, 
pH 7.5, 50 MM DHQ 40 |iM compound, 1 nM «fl3HQD). 

Inhibitor IC 50 and Mode of Inhibition Experiments 

Inhibitor analogs were purchased from Chembridge Corpora- 
tion and maintained as 50 mM stocks in DMSO. IC 5 o measure- 
ments were performed in triplicate using the primary screen assay 
and varying concentrations of inhibitor (0-400 |lM). Reaction 
rates were monitored continuously for up to 30 minutes, and 
slopes were calculated from the linear portions of the progress 
curves and corrected for background signal. Data were fit to 
equation 1 

v = a /(l + [7]/7C 50 ) (1) 

using the Enzyme Kinetics module of SigmaPlot (v. 9.01 Systat 
Software, Inc.) where v is the reaction rate in the presence of 
inhibitor, a is the reaction rate in the absence of inhibitor, and [I] 
is the inhibitor concentration. Mode of inhibition studies with 
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Figure 1. The high-throughput screen assay used to identify the three cdDHQD inhibitors. (A). The conversion of 3-dehydroquinate 
(DHQ) to 3-dehydroshikimate (DHS) by cdDHQD is coupled to saturating concentrations of two downstream enzymes in the shikimate pathway, 
shikimate dehydrogenase (SDH) and shikimate kinase (SK). The loss of NADPH fluorescence in the SDH-catalyzed step of the reaction provides a 
measurable, continuous readout for the assay. The presence of SK drives the reversible SDH reaction to completion and increases the dynamic range 
of the assay. (B) The structure of the three hit compounds and corresponding K, values. 
doi:1 0.1 371 /journal.pone.0089356.g001 



compounds 1 and 3 were performed in duplicate using the 
primary screen assay with varying concentrations of inhibitor and 
substrate. Background measurements of the assay components 
lacking «/DHQD were performed in duplicate at all substrate and 
inhibitor concentrations. Reaction rates were monitored contin- 
uously for up to 30 minutes, and slopes were calculated from the 
linear portions of the progress curves and corrected for 
background signal. The data were fit to the equations describing 
competitive, non-competitive, mixed, and uncompetitive inhibi- 
tion models using the Enzyme Kinetics module of SigmaPlot. 

Fluorescence Thermal Shift (FTS) Assays 

To establish dose-response relationships, compounds 1-3 were 
assayed in quadruplicate at concentrations ranging from 2.5- 
320 mJVL Compounds 1-3 were added to 1.5 [Ig of i/DHQD, 
«/DHQD, or «DHQD and 10 nL of 5000x Sypro Orange 
(Invitrogen) in 10 (XL of assay buffer (10 mM Tris [pH 8.4] and 
150 mM NaCl). After compound addition, the plate was sealed 
with an optical seal, shaken, and centrifuged. The CFX384 (Bio- 
Rad Laboratories) real-time PCR machine was used for fluores- 
cence detection and thermal scanning, which progressed from 10 
to 95°C with a temperature ramp rate of 1.5°C/min. For analysis, 
the midpoint of thermal denaturation (T m ) for each condition was 
calculated using an in-house software package. 

NMR Methods 

NMR experiments were performed on a Bruker 900 MHz 
AVANCE spectrometer equipped with a cryogenic triple reso- 
nance probe. Experimental conditions were 10-15 |jM DHQD, 
10-500 m_M compound in PBS/pH 7.2 in 90% 'H 2 0, 10% 2 H 2 0 
(WaterLOGSY) or 100% 2 H z O (STD) at 25°C. The Water- 
LOGSY experiments were performed as previously described [1 1— 



1 3] . Water was selectively saturated using a 2 msec square shaped 
pulse with a mixing time of 1-2 sec and a relaxation delay of 
2.5 sec. STD experiments were performed as previously described 
[1 1,14, 15]. In these experiments, protein 'H were saturated with a 
train of 50 msec gaussian-shaped pulses at 100 Hz power for 1 sec 
with "on" resonance saturation at — 1 ppm and "off' resonance 
saturation at 30 ppm (the relaxation delay was 2.5 sec before the 
saturating pulses). Spectra were processed by NMRPipe with a 
5 Hz line broadening function and analyzed by NMRDraw [16]. 
Relative %STD was defined as 100 X STD obs /STD milx where 
STD = AI/I off , AI = I off — I on and l 0 g and I on are the intensities 
observed for the various resonances after the "off' and "on" 
presaturation of HA. Errors in the WaterLOGSY and STD were 
estimated as AI/I re£ ((NAi/AI) 2 +(N Iofr /I rrf ) 2 ) () ' 5 [14], where NAI and 
N Iref are the noise calculated by NMRDraw in the appropriate 
spectrum (no protein in WaterLOGSY and I off in the STD). 
Titration data were fit to the appropriate functions using 
Kaleidagraph 4.1.3. 

Results 

Discovery of cdDHQD Inhibitors 

The dehydration of DHQ, to DHS can be monitored 
spectroscopically at 234 nm due to the increased absorption of 
the product [17]. However, due to signal interference from library 
compounds at this wavelength such an assay is poorly suited for 
HTS. Therefore, we coupled the DHQD reaction to shikimate 
dehydrogenase, the subsequent step of the shikimate pathway and 
monitored the conversion of NADPH to NADP using a 
fluorescence-based assay (Figure 1A). During assay development 
we noted that the reaction was not progressing to completion and, 
therefore, shikimate kinase, the enzyme that follows shikimate 
dehydrogenase in the shikimate pathway, was added to shift the 
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equilibrium towards the reaction product. This assay design gave 
an improved signal and Z'-factor of 0.68, a value indicative of a 
sensitive and robust assay [18], 

To identify small molecule inhibitors of cdDHQD, a high 
throughput screen (HTS) of 50,000 Chembridge compounds was 
performed. This assay generated a low hit rate, with the vast 
majority of false positives resulting from compound interference 
with the fluorescent signal and a minority of false positives 
affecting the coupling enzymes. Using the direct UV assay to 
validate hits identified by HTS, we discovered three cdDHQD 
inhibitors (compounds 1, 2, and 3) with ~30 uM IC 50 values 
(Figure IB). 

Hit compounds 1 and 2 exhibit striking structural similarity, 
both being 2,5-dichlorobenzenesulfonamide derivatives. While 
compound 1 incorporates a 6-substituted 2-aminobenzo[d]thia- 
zole group, compound 2 has a 5-phenylthiazol-2-amine substitu- 
ent. Interestingly, compound 3, an amide of tetrahydroben- 
zo[A]fhiophene-3-carboxylic acid, is structurally unrelated to 1 and 
2. All compounds appear to provide good starting points for future 
optimization in the hit-to-lead process. 

Kinetic Characterization of Hit Compounds 

To better understand modes of inhibition and to determine 
inhibition constants (Kj) for the three hits, we performed a 
complete analysis using steady state kinetics. Initial reaction rates 
measured at different substrate (DHQ) and compound 1 
concentrations can be fit using either competitive or mixed-type 
inhibition models (Figure S1A in File SI). The K ; value of 
compound 1 to the cdDHQD was calculated to be 11 to 15 uM, 
depending on the kinetic model used. Low solubility of compound 
2 prevented us from accurately determining its K;. Ambiguity 
between kinetic models was also observed for compound 3, for 
which K ; values ranged between 19 and 22 |iM depending of the 
kinetic model used (Figure SIB in File SI). 

In the competitive model, the inhibitor (I) binds to the enzyme 
(E) but not the enzyme-substrate (ES) complex - only EI is formed. 
In contrast, in the mixed inhibition model, the inhibitor binds to 
both E and ES complex, forming EI and ESI. Competitive 
inhibition is characterized by an increased apparent K m but 
unchanged V max , as increased substrate concentration will shift EI 
to ES, negating the inhibition. For mixed inhibition, while the 
effect due to the presence of the inhibitor on apparent K m can 
vary, V max is decreased - even at high [S] , the presence of ESI 
remains, reducing V max . Thus, this differential effect on V max can 
be used to distinguish between competitive and mixed models of 
inhibition. To clarify the mechanism of inhibition, we re-assayed 
cdDHQD with substrate concentrations reaching as high as 5 mM, 
asking if the V max differs to that observed in the absence of 
inhibitor. Indeed, in the presence of 200 U.M of compound 3 the 
apparent K m increased 10-fold but V max was unchanged (results 
not shown). This result is consistent with a competitive mechanism 
of inhibition and is counter to the mixed inhibition model. 

Selectivity of Compounds 1-3 

To determine inhibitor selectivity, we screened compounds 1-3 
against a panel of purified DHQDs that included the type II 
btDHQD, type II acDHQD, type II j/DHQD, and type I 
.rfDHQD. &DHQD shares 56% sequence identity (78% sequence 
homology) and a highly homologous quaternary structure 
(backbone root mean square deviation of homodimer= 1.2 A) 
with cdDHQD [19]. As anticipated, the kinetic assays revealed that 
the compounds were inactive against the structurally and 
mechanistically distinct type II DHQDs, but, surprisingly, also 
inactive against the related type I feDHQD (Table 1). 



Table 1. 


IC 50 s for Multiple DHQDs (uM). 








Compound 1 


Compound 2 


Compound 3 


Type 1 


cdDHQD 


34 


35 


31 


seDHQD 


>200 


>200 


>200 


Type II 


btDHQD 


>200 


>200 


>200 


vcDHQD 


>200 


>200 


>200 


ypDHQD 


>200 


>200 


>200 



doi:1 0.1 371 /journal.pone.0089356.t001 



Next, we turned to two biophysical techniques to confirm the 
unexpected intra-type selectivity. First, Fluorescence Thermal Shift 
(FTS) assays were performed to probe the effect of compounds 1 
3 on caDHQP, btDHQD, and kDHQD thermal stability. 
Consistent with selective cdDHQD binding, all three compounds 
thermally stabilized cdDHQD, but not btDHQD and .seDHQD, in 
a dose-dependent manner (Figure 2). Second, NMR spectroscopy 
was used to measure compounds 1-3 resonances in the presence 
and absence of enzyme. In the case of compounds 1 and 2, we 
observed line broadening of the compound resonances upon 
addition of caDHQD, which is characteristic of chemical exchange 
in the intermediate time scale (Figure S2A in File SI). The 
double reciprocal plot of the line broadening effects yields a K d of 
—27 uM for binding of compound 1 and of ~25 |XM for 
binding of compound 2 to cdDHQD (Figure 3A & 3B). 
Importantly, little line broadening was observed for reDHQD at 
similar compound concentrations (<1 Hz vs. >6 Hz for 
«DHQD, Figure S2B), preventing the calculation of a Kj and 
suggesting that compounds 1 and 2 bind with significantly lower 
affinity to reDHQD. 

Due to a faster exchange rate, the resonances of compound 3 
were not measurably broadened upon the addition of cdDHQD. 
Accordingly, we chose the WaterLOGSY NMR experiment, 
which has previously been shown to be useful for small molecule 
binding to proteins in fast exchange [11-13], to assay binding of 
compound 3 to DHQD. The dependence of the WaterLOGSY 
signal on concentration allowed us to determine a Kh of ~65 U.M 
for binding of compound 3 to cdDHQD (Figure 3C). Impor- 
tantly, in the case of teDHQD, the affinity for compound 3 is again 
significandy lower (Figure 3D and Figure S3 in File SI), 
confirming its selectivity for cdDHQD. 

NMR Competition Experiments Suggest a Shared Binding 
Site for Compounds 1-3 

Due to the structural similarity between compounds 1-3, we 
next tested whether they bind to overlapping sites on cdDHQD 
using an NMR-based competition assay. This experiment probes 
how inhibitors affect product binding to the enzyme. We first 
incubated the enzyme with substrate allowing the reaction to reach 
equilibrium - the reaction is monitored by comparing the relative 
intensity of substrate and product resonances. In this way we 
determined that at equilibrium, the ratio between substrate/ 
product to be ~ 1 :20, a result consistent with literature reports that 
also found a high equilibrium ratio [20]. Since the substrate to 
product ratio is skewed towards the products at equilibrium 
conditions we decided to evaluate the effect of inhibitors on 
product resonances. 
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Figure 2. Fluorescence Thermal Shift demonstrates that Compounds 1-3 bind selectively to caDHQD. (A) The effect of 160 uM DHS/ 
compounds 1-3 on the thermal stability of cdDHQD, seDHQD, and bfDHQD as determined by FTS. The AT m represents the difference in midpoint 
thermal denaturation temperature compared to the ligand-less control. (B) FTS experiments show AT m of cdDHQD as a function of compound 1-3 
concentration. 

doi:1 0.1 371 /journal.pone.0089356.g002 



Hence, prior to addition of the compounds, we incubated the 
enzyme in the presence of substrate until equilibrium was 
established, and subsequently monitored the resonances of the 
product DHS. As shown in Figure 4, the addition of compounds 
1-3 reduces DHS binding by 40-50%. This result is consistent 
with all three compounds binding to the product binding site on 
tt/DHQD (primary data are presented in Figure S4 in File SI). 



However, this experiment by itself does not exclude the possibility 
of the inhibitors binding to an allosteric binding site, potentially 
even different allosteric sites for each compound. 




0 50 100 150 200 0 100 200 300 400 500 

compound 3 uJVl compound 3 u.M 



Figure 3. NMR demonstrates that Compounds 1-3 bind selectively to cdDHQD. NMR titration of (A) compound 1, (B) compound 2, (C) 
compound 3 to cdDHQD. All three compounds demonstrate binding to cdDHQD with Kd values of —25 U.M for compounds 1 & 2, and —65 U.M for 
compound 3. (D) NMR titration of compound 3 to seDHQD reveals much reduced affinity, with a K d of —400 |iM binding. In (A) and (B), the red and 
black lines represent the results from two different resonances. 
doi:1 0.1 371 /journal.pone.0089356.g003 
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Sample 

Figure 4. NMR competition experiments to characterize the 
binding site of compounds 1-3 on cdDHQD. WaterLOGSY NMR 
competition experiment between compounds 1-3 and product DHS in 
the presence of cc/DHQD. The reduced WaterLOGSY signal for DHS in 
the presence of the compounds indicate that the compounds compete 
with DHS for binding to cc/DHQD. 
doi:1 0.1 371 /journal.pone.0089356.g004 

NMR Experiments to Map the Compound-DHQD 
Interaction 

Saturation Transfer Difference (STD) NMR has been shown to 
be useful for identifying the ligand 'H in closest proximity to the 
protein surface [11,14,15], Accordingly, we used the technique to 
characterize the interactions of compounds 2 and 3 with «/DHQD 
(compound 1 was not amenable to this assay due to line- 
broadening and weak signal). A summary of these experiments is 
shown in Figure 5. In the case of compound 2, the relative STD 
signals range from 57-100%, with protons from the dichloroben- 
zenesulfonamide substituent producing reliably stronger signals 
then those from the phenylthiazol core. In the case of compound 
3, the relative STD signals range from 41-100%, with the 
thiophene proton producing the strongest signal, followed by 
protons on the phenyl ring. Interestingly, protons on compound 
2's dichlorobenzenesulfonamide substituent and compound 3's 
phenyl group, as well as compound 2's phenylthiazol and 
compound 3's tetrahydrobenzo [b] thiophene group, exhibit com- 
parable relative intensities, raising the possibility that the two 
compounds similarly engage cdDHQT). 

Structure-activity Analysis 

To identify the specific elements required for DHQD inhibition, 
we tested analogs of compounds 1-3 (Figure S5A in File SI). 

For compounds 1 and 2, modifications to the dichlorobenzene- 
sulfonamide core were poorly tolerated. Analogs lacking the two 
chlorines were completely inactive and replacing the chlorines 
with fluorines resulted in much reduced inhibition. Adding a third 




chlorine at the ortho position also abolished inhibition suggests 
that the dichlorobenzenesulfonamide core makes critical and easily 
disrupted interactions with the enzyme. As is evident from the 
differences between compounds 1 and 2, modifications to the 
thiazole substituent were better tolerated. While a compound 
containing the unmodified thiazole ring was inactive, a variety of 
compounds with distinct additions to the thiazole substituent 
displayed similar activity to the preliminary hits (Figure S5B in 
File SI) 

Discussion 

To our knowledge, the compounds described here represent the 
first non-mechanism based type I DHQD inhibitors identified to 
date. Moreover they exhibit properties of inter- and intra-type 
specificity desirable of a narrow-spectrum antibiotic for CDI 
treatment. While the identification of selective «/DHQD inhibitors 
provides a proof-of-concept, additional studies are required to 
optimize inhibitor binding. In principle, the present set of studies 
provides insight that should be useful for the development of high 
affinity inhibitors. In particular, the STD-NMR data suggests that 
compound 1 and 2's dichlorophenyl core intimately associates 
with crfDHQD and SAR deduced from the analysis of inhibitor 
analogs confirm its importance for binding. By contrast, these 
approaches suggest that modified thiazole substituent is less 
intimately associated with the protein and less critical for binding. 
These observations are consistent with the existence of a region 
adjacent to the thiazole binding site that can be viably expanded 
into - arguably making substitutions at the thiazole ring to be the 
most promising strategy for inhibitor optimization. 

It is clear the structural information about inhibitor binding 
could help explain the observed specificity and further aid the 
development of high affinity inhibitors. Despite attempting 
multiple co-crystallization and soaking experiments, we were 
unable to obtain inhibitor bound crystal complexes. In the absence 
of a complex crystal structure, inferences about the basis of 
inhibitor binding can be made from available structural and 
kinetic data. Previously described crystal structures reveal DHQD 
adopts two discrete conformational states [21]. In the absence of 
ligand the functionally important (38-a8 loop adopts an open and 
partially disordered conformation. In liganded structures, closure 
of the (38-0(8 loop establishes hydrogen bonding interactions with 
the lysinel71 Schiff base-bound reaction intermediate [19]. 

Three observations suggest that the inhibitors, unlike the 
reaction intermediate, bind DHQD's open loop conformational 
state: First, whereas co-crystallization and crystal soaking exper- 
iments readily yielded reaction intermediate bound complexes, to 
date, we have been unable to obtain inhibitor bound complexes, as 
noted above (data not shown). If the substrate and inhibitor bind 
the same conformational state, then successful soaks with the 
comparable affinity inhibitors could reasonably be expected. 




Figure 5. STD NMR experiments map the binding of compounds 2 and 3 on cdDHQD. The relative proximity of compounds' carbons to 
cc/DHQD atoms based on STD NMR data is shown. 
doi:1 0.1 371 /journal.pone.0089356.g005 
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Second, with loop closure the reaction intermediate is buried 
beneath the protein surface [2 1] . As the active site snuggly fits the 
reaction intermediate, it is difficult to envision how it could 
accommodate the appreciably larger structures of compounds 1— 
3. Finally, eaDHQD and reDHQD active sites are virtually 
structurally identical in the closed (38-a8 loop conformational 
state. Such structurally similar binding sites are inconsistent with 
the drastically differing binding affinities the compounds exhibit 
for the two enzymes. 

Collectively, these observations make a strong case that 
compounds 1—3 do not bind the closed loop conformational state. 
The size of the compounds necessitates that they bind the open 
loop conformational state and extend out of the active site. This 
fact could account for the specificity of the compounds to 
caDHQD. A number of residues on the perimeter of the active 
site differ between eaDHQD and nDHQD and interactions with 
one or more of these likely accounts for the observed intra-type 
specificity. Future work will focus on testing the ability of these 
compounds to inhibit the growth of C. difficile, while at the same 
time, sparring a large subset of the commensal flora. 

In sum, this work demonstrates that selective inhibition of 
caDHQD is possible, and that the hits identified bind at the active 
site, and in so doing, likely prevent the (38-oc8 loop from closing. 
The interactions made with residues more distant to the active site 
are likely responsible for the selectivity towards the C. difficile 
enzyme versus the highly homologous S. enterica enzyme. We are 
currently synthesizing analogs of compounds 1-3 in order to 
identify molecules that bind with higher affinity, but that also 
maintain a preference for cdDHQD. Such compounds would hold 
promise for CDI therapy. 

Supporting Information 

File SI This file contains Figure Si-Figure S5. Figure SI, 
Fit of initial rate measured under varying substrate and inhibitor 
concentrations to the mixed-type and competitive models. (A) Fits 
for compound 1. (B) Fits for compound 3. Figure S2, NMR data 

References 

1 . Rupnik M, Wilcox MH, Gerding DN (2009) Clostridium difficile infection: new 
developments in epidemiology and pathogenesis. Nat Rev Microbiol 7: 526-536. 

2. Kuijper EJ, van Dissel JT, Wilcox MH (2007) Clostridium difficile: changing 
epidemiology and new treatment options. Curr Opin Infect Dis 20: 376-383. 

3. Poutancn SM, Simor AE (2004) Clostridium difficile-associated diarrhea in 
adults. CMAJ 171: 51-58. 

4. Mullanc KM, Miller MA, Weiss K, Lcntnck A, Golan Y, ct al. (201 1) Efficacy of 
lidaxomicin versus vancomycin as therapy for Clostridium difficile infection in 
individuals taking concomitant antibiotics for other concurrent infections. Clin 
Infect Dis 53: 440-447. 

5. Crcdito KL, Appelbaum PC (2004) Activity of OPT-80, a novel macrocycle, 
compared with those of eight other agents against selected anaerobic species. 
Antimicrob Agents Chemother 48: 4430-4434. 

6. Herrmann KM (1995) The shikimate pathway as an entry to aromatic secondary 
metabolism. Plant Physiol 107: 7-12. 

7. Gourley DG, Shrive AK, Polikarpov I, Krell T, Coggins JR, ct al. (1999) The 
two types of 3-dehydroquinase have distinct structures but catalyze the same 
overall reaction. Nat Struct Biol 6: 521-525. 

8. Klcanthous C, Deka R, Davis K, Kelly SM, Cooper A, ct al. (1992) A 
comparison of the enzymological and biophysical properties of two distinct 
classes of dehydroquinase enzymes. BiochemJ 282 (Pt 3): 687—695. 

9. Butler J, Alworth W, Nugent M (1974) Mechanism of dehydroquinase catalyzed 
dehydration. I. Formation of a Schiff base intermediate. JAmChemSoc 96: 
1617-1618. 

10. Harris JM, Gonzalcz-Bcllo C, Klcanthous C, Hawkins AR, Coggins JR, ct al. 
(1996) Evidence from kinetic isotope studies for an enolate intermediate in the 
mechanism of type II dehydroquinases. BiochemJ 319 (Pt 2): 333-336. 

1 1 . Antanasijcvic A, Cheng H, Wardrop D, Rong L, Caffrcy M (20 1 3) Inhibition of 
Influenza H7 Hemagglutinin-Mediated Entry. PLoS One: in press. 



for binding of compounds 1 and 2 to «/DHQD. (A) Line 
broadening observed for compound 1 upon addition of osDHQD 
or seDHQD. (B) Line broadening observed for compound 2 upon 
addition of eaDHQD or yeDHQD. Note that line broadening was 
observed only in the presence of cdDYLQT) but not with .seDHQD, 
demonstrating the selectivity of the compounds to o/DHQD. The 
experimental conditions were 10 |lM DHQD, 100 |lM compound 
in PBS/pH 7.2, 10% 2 H 2 0 at 25°C. Data were acquired on a 
Bruker 900 MHz spectrometer equipped with a cryogenic triple 
resonance probe. Line widths were estimated using NMRDraw 
(Delaglio et al., 1995). Figure S3, NMR data for binding of 
compound 3 to caDHQD. The experimental conditions were 
10 MM DHQD, 100 |JM compound 3 in PBS/pH 7.2, 10% 
2 H 2 0 at 25°C. Data were acquired on a Bruker 900 MHz 
spectrometer equipped with a cryogenic triple resonance probe. 
The WaterLOGSY experiments employed a 1 sec mixing time. 
Figure S4, NMR data for the WaterLOGSY competition assay. 
Experimental conditions were 10 |JM DHQD, 100 |lM DHS and 
100 |0.M compounds 1, 2 or 3 in PBS/pH 7.2, 10% 2 H 2 0 at 
25°C. Data were acquired on a Bruker 900 MHz spectrometer 
equipped with a cryogenic triple resonance probe with a 
WaterLOGSY mixing time of 2 sec. Figure S5, Structure-activity 
relationship analysis for compounds 1—3. (A), (B), and (C): 
Inhibition of the a/DHQD catalyzed dehydration reaction was 
measured for compounds 1-3 and select analogs. The critical role 
for binding through the phenyl aryl chlorine atoms (compounds 1 
& 2) and the carboxamideylic moiety (compound 3) is apparent. 
(D) A series of analogs demonstrate that distinct additions to the 
dichlorobenzenesulfonamide thiazole core shared by compounds 1 
and 2 have little effect on IC 50 . 
(DOCX) 

Author Contributions 

Conceived and designed the experiments: AL KR SL AA WA MC. 
Performed the experiments: KR SHL AA. Analyzed the data: KR SHL 
AA MC AL. Wrote the paper: SHL AL. 



12. Dalvit C, Fogliatto G, Stewart A, Vcroncsi M, Stockman B (2001) Water- 
LOGSY as a method for primary NMR screening: practical aspects and range of 
applicability. J Biomol NMR 21: 349-359. 

13. Ramirez B, Antanasijcvic A, Caffrcy M (2013) Ligand screening using NMR. 
Meth Molecul Biol: in press. 

14. McCullough C, Wang M, Rong L, Caffrey M (2012) Characterization of 
influenza hemagglutinin interactions with receptor by NMR. PLoS One 7: 
c33958. 

15. Meyer B, Peters T (2003) NMR spectroscopy techniques for screening and 
identifying ligand binding to protein receptors. Angew Chcm Int Ed Engl 42: 
864-890. 

16. Delaglio E, Grzcsick S, Vuister GW, Zhu G, Pfeifer J, et al. (1995) NMRPipe: a 
multidimensional spectral processing system based on UNIX pipes. J Biomol 
NMR 6: 277-293. 

17. Chaudhuri S, Lambert JM, McColl LA, Coggins JR (1986) Purification and 
characterization of 3-dehydroquinase from Escherichia coli. Biochem J 239: 
699-704. 

18. ZhangJH, Chung TD, Oldenburg KR (1999) A Simple Statistical Parameter for 
Use in Evaluation and Validation of High Throughput Screening Assays. 
J Biomol Screen 4: 67—73. 

19. Light SH, Minasov G, Shuvalova L, Duban ME, Caffrcy M, ct al. (2011) 
Insights into the mechanism of type I dehydroquinate dehydratases from 
structures of reaction intermediates. J Biol Chcm 286: 3531-3539. 

20. Mitsuhashi S, Davis BD (1954) Aromatic biosynthesis. XII. Conversion of 5- 
dchydroquinic acid to 5-dehydroshikimic acid dy 5-dehydroquinase. Biochim 
Biophys Acta 15: 54—61. 

21. Light SH, Minasov G, Shuvalova L, Peterson SN, Caffrcy M, ct al. (2011) A 
conserved surface loop in type I dehydroquinate dehydratases positions an active 
site arginine and functions in substrate binding. Biochemistry 50: 2357-2363. 



PLOS ONE | www.plosone.org 



7 



February 2014 | Volume 9 | Issue 2 | e89356 



